function out = T2_fun3(beta_d,lambda_c,muv,I)
% muv must be row vector I-dim row vector

temp0 = lambda_c+(I+1)*beta_d;

out2 = -(temp0*muv - beta_d*sum(muv))/(temp0*(beta_d+lambda_c)) + muv/(lambda_c+2*beta_d);
out1 = -beta_d*sum(out2);

out = [out1 out2];  % output is (I+1)-dim row vector
